Dynamics-dependent criticality in models with q absorbing states 
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We study a one-dimensional, nonequilibrium Potts-like model which has q symmetric absorbing 
states. For g = 2, as expected, the model belongs to the parity conserving universality class. For 
5 = 3 the critical behaviour depends on the dynamics of the model. Under a certain dynamics 
it remains generically in the active phase, which is also the feature of some other models with 
three absorbing states. However, a modified dynamics induces a parity conserving phase transition. 
Relations with branching-annihilating random walk models are discussed in order to explain such a 
behaviour. 
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I. INTRODUCTION 



Recently, dynamical and nonequilibrium properties of 
many-body systems have been intensively studied. Of 
particular interest are nonequilibrium phase transitions 
which might appear in the stationary state of such sys- 
tems [0. It is believed that continuous phase transitions 
can be classified into relatively few universality classes. 
It is becoming evident, however, that such a classifica- 
tion is far more complicated in non-equilibrium systems 
than in equilibrium ones. In the equilibrium case, renor- 
malization group and conformal-invariance theories gave 
a powerful description of this universal phenomena. In 
non-equilibrium, the situation is far less clear due to the 
lack of a general theory. Accordingly, a large body of 
works is based on numerical simulations. 

Models with absorbing states constitute a particularly 
rich class. For these models there is a substantial ev- 
idence that continuous phase transitions can be classi- 
fied into some universality classes. In particular, a large 
group of models falls into the so-called directed percola- 
tion universality class (DP) and it was conjectured that 
all models with a single absorbing state, positive one- 
component order parameter and short-range dynamics 
should generically belong to this universality class 
Another group consists of models having a double de- 
generate absorbing state or whose dynamics obey some 
conservation law; they belong to the parity-conserving 
universality class (PC) 

The behaviour of models with a larger number of ab- 
sorbing states was also addressed in the literature. For 
example Bassler and Browne examined a model with 
three absorbing states and concluded that depending on 
some parameters their model might exhibit DP or PC 
criticality However, the control parameters in their 
model (adsorption rates) introduce asymmetry between 
absorbing states when the critical point is approached. 
Consequently, one or two species are effectively expelled 
from the system upon approaching the critical point 



hence the DP or PC criticality is an expected feature 
of this model. Similarly, certain asymmetries are respon- 
sible for the DP criticality in yet another model with mul- 
tiple absorbing states which was studied by Janssen |5| . 

To examine the role of degeneracy, one has to study 
models where the symmetry between absorbing states is 
not broken at the level of dynamics. Good candidates for 
such a system are certain multi-species generalizations 
of the contact process model In the two-species 

case these models exhibit the expected PC criticality. In 
the three-species case Hinrichsen suggested Q that such 
models will always remain in the active phase. This prop- 
erty, which should also be true for models with a larger 
number of absorbing states, follows from an approxi- 
mate relationship with the g-species, parity conserving 
branching- annihilating random walk models (g-BARW2, 
where the number at the end indicates the number of 
offsprings) Recently, Hinrichsen's conjecture was nu- 
merically confirmed for models with three and four ab- 
sorbing states 1^. 

One can thus expect that for a given dimensional- 
ity, the number of (symmetric) absorbing states is the 
relevant parameter determining the critical behavior of 
a given model. However, for nonequilibrium systems 
some other details of the dynamics like e.g., exclusion 
effects or certain local symmetries might 

affect critical behaviour. The goal of the present pa- 
per is to provide yet another example of the dynamics- 
dependent criticality. In particular, we show that for 
a one-dimensional model with three absorbing states its 
critical behavior is PC-like (instead of the expected q- 
BARW2-like), when certain local constraints on the dy- 
namics of a model are introduced. Let us emphasize that 
these constraints, which can be regarded as a local sym- 
metry breaking, do not violate the symmetry between ab- 
sorbing states [0. We also suggest a mechanism which 
could explain such a behaviour. 

Our work show that not only global properties of a 
given model, like the number of absorbing states, but 
also certain details of the dynamics are relevant to deter- 
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mine its critical behaviour. Accordingly, classification of 
the critical properties of models with absorbing states is 
more complicated than originally thought. 

The critical behaviour of g-BARW2 models was re- 
cently found to be richer than originally expected. In- 
deed, in one dimension hard-core effects are known to 
change the off-critical exponents of the model ||ll|,|l^ . Let 
us notice that these effects do not change the location of 
the critical point and of the on-critical exponents . At 
the coarse-grained level, our model can be also regarded 
as a certain g-BARW2 model. In our case, however, dy- 
namical details have more dramatic effect: they change 
both the location of the critical point and values of all 
critical exponents. 

II. MODEL AND ITS MONTE CARLO 
SIMULATIONS 

Before defining our dynamical model, let us recall some 
basic properties of the usual equilibrium Potts model. 
First we assign at each lattice site i a g-state variable 
CTi = 0, 1, ...,q — 1. Next, we define the energy of this 
model through the Hamiltonian: 

H = -Y,S,^,^, (1) 

where summation is over pairs of i and j which are usu- 
ally nearest neighbours and S is the Kronecker delta func- 
tion. This equilibrium model was studied using many 
different analytical and numerical methods and is a rich 
source of the information about phase transitions and 
critical phenomena 

To simulate numerically the equilibrium Potts model 
defined using the Hamiltonian (|l|) one introduces a 
stochastic Markov process with transition rates chosen in 
such a way that the asymptotic probability distribution 
is the Boltzmann distribution. One possibility of choos- 
ing such rates is the so-called Metropolis algorithm. In 
this method [|6j one looks at the energy difference AE 
between the final and initial configuration and accept the 
move with probability min {1,6-^^^/^}, where T is tem- 
perature measured in units of the interaction constant of 
the Hamiltonian (|^), which was set to unity. To obtain 
a final configuration one selects randomly a site and its 
state (one out of q in our case). In the above described 
algorithm for T > there is always a positive probability 
of leaving any given configuration (even when the final 
configuration has a higher energy). Accordingly, such a 
model does not have absorbing states for T > 0. 

A. A-model 

To transform the standard Metropolis dynamics into 
the dynamics with absorbing-states we make the follow- 



ing modification: 
Restriction A: 

When all neighbours of a given site are in the same state 
as this site, then this site cannot change its state (at least 
until one of its neighbours is changed) . 

In the following, this nonequilibrium model will be re- 
ferred to as A-model. Obviously, any of q ground states 
of model (^ is an absorbing state of A-model and the 
dynamics does not favour any of the absorbing states. 
Since the dynamics of our models is obtained from a 
modification of the Metropolis algorithm of an equilib- 
rium system, transition probabilities are parametrized 
by temperature-like quantity T. Strictly speaking, for 
our model the ordinary (i.e., equilibrium) temperature 
cannot be defined. Nevertheless, we will refer to this 
quantity as temperature. 

To study the properties of this A-model we performed 
Monte Carlo simulations. A natural characteristic of 
models with absorbing states is the steady-state density 
of active sites p. A given site i is active when at least one 
of its neighbours is in a state different than i. Otherwise 
the site i is called nonactive. In addition to the steady- 
state density we also looked at its time dependence p(t). 
In the active phase p{t) converges to the positive value 
while at criticality, p{t) has a power-law decay p ^ t~^ . 
In the absorbing phase the density p decays either faster 
than a power of t or as power of t but with a different 
exponent than the critical exponent 5. The unit of time 
is defined as a single (on average) update of each active 
site. 

In addition, we used the so-called dynamic Monte 
Carlo method where one sets the system in the absorb- 
ing state, locally initiate activity, and then monitor some 
stochastic properties of surviving runs . The most fre- 
quently used characteristics are the survival probability 
P{t) that the activity survives at least until time t and 
the number of active sites N{t) (averaged over all runs). 
At criticality these quantities are expected to have power- 
law decay: P{t) - t~^' and N{t) - VK (For some models 
5 = 5', but exceptions are also known [1|). 

Our simulations were made for various system sizes and 
we ensured that the system was large enough so that our 
results are size independent. 

The simplest case is to consider a one-dimensional 
chain. However, for q = 2 and for any temperature T, 
this model is trivially equivalent to the T = tempera- 
ture Ising model with Metropolis dynamics. Indeed, in 
this case the allowed moves are only those which do not 
increase energy and they are always accepted. The same 
rule governs the dynamics of the T = Ising chain. 

To overcome this geometrical 'pathology' we consider 
our A-model on a ladder-like lattice where two chains 
are connected by inter-chain bonds such that each site 
has three neighbours. Fig. |^ illustrates the possibil- 
ity of spreading of activity in the ladder geometry for 
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q = 2. The results of the simulations of such case are 
in agreement with the up-to-date knowledge concern- 
ing one-dimensional systems with multi-absorbing states. 
Thus we only briefly describe our results. For q = 2 the 
model has two symmetric absorbing states (all ai — or 
1). Qualitatively, in this case the model resembles other 
models with two absorbing states , which are known 
to belong to PC universality class and our simulations 
confirm such a behaviour. 



X 


X 


X 











X 


X 


X 











X 


X 


X 


X 








X 


X 


X 











X 


X 


X 


X 








X 


X 














X 


X 


X 


X 








X 


X 





X 








X 


X 





X 








X 


X 





X 









FIG. 1. At positive temperature a single domain can 
branch into additional domains. Successive (in time) con- 
figurations differ only by a single-site flips. States and 1 are 
represented by "x" and "o", respectively. In the single-chain 
geometry and q = 2 such a branching is forbidden and domain 
walls can only diffuse or annihilate. In terms of BARW mod- 
els, where each domain wall represents a particle, the above 
sequence is equivalent to branching. 



active sites suggests a continuous phase transition be- 
tween active and absorbing phases around T = 2.7 (see 
Fig. 1^). Measuring the time dependence of p{t) wc ob- 
served that at T = 2.7 p{t) ~ where S is very close to 
the PC value 0.286. Moreover, in the low-temperature 
(T < 2.7) phase p{t) ~ t^"-^, which is also a typi- 
cal feature of PC models. Additional confirmation of 
PC criticality in this case is obtained using the dynam- 
ical Monte Carlo method which yields S' = 0.29(3) and 
r] = -0.02(3). 

A different behaviour appears in the q = 3 case. Al- 
though in Fig. H one can see a sudden change of the order 
parameter around T = 0.8, there is no phase transition 
in this case. Examining the behaviour of p{t), we checked 
that even at low temperature (T = 0.5 and 0.6) the sys- 
tem remains in the active phase. In addition simulations 
suggest that for q — 4: and 5 A-model behaves similarly 
to the q — 3 case. 

Actually, we expect that for q > 3 the model has a 
critical point but only at T = 0. This point corresponds 
to the case of zero branching rate in a (3'-BARW2 model, 
which is known to be characterized by e.g., /3 = 1 (with 
the order parameter expressed in terms of the branching 
rate) §. 

As we already mentioned, the absence of the transition 
for 9 > 3 and T > is an expected feature. However, 
as we will show below, additional restriction in the dy- 
namical rules of our model induces a transition even for 
(J > 3. 



B. B-model 
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FIG. 2. The steady-state density p as a function of T 
for A-model model and q = 2 {+), 3 (*) and for B-model 
with q — 3 (□). Simulations were made for the system 
size L = 30000 and simulation time t = W^MCS. The 
low-temperature tail for q = 3 does not diminish for longer 
simulations or larger system size. 

First, the steady-state measurement of the density of 



This restriction can be formulated as follows: 
Restriction B: 

A flip into a state different than any of its neighbours is 
prohibited. In other words, spontaneous creation of for 
example domains of type A between domains of type B 
and C is forbidden. Here, A, B, and C denote three (out 
of q) different states. 

Let us notice that restriction B satisfies at the same 
time restriction A (but of course not vice versa). In the 
following we will refer to the model satisfying restric- 
tion B as B-model. Let us also notice that restriction B 
does not break the symmetry and B-model similarly to 
A-model has q symmetric absorbing states. Results of 
the simulations of B-model for g = 3 are shown in Fig. ^ 
Fig. ^ (for q = 2 both dynamics A and B are equiva- 
lent). First, let us notice that the density p (Fig. ||) is 
only slightly larger for B-model than for A-model with 
q = 2 (the difference is, however, larger than error bars). 
In addition, this difference diminishes upon approaching 
the transition point and within our numerical accuracy 
both transitions seem to take place at the same temper- 
ature. Later on we provide some arguments which could 
explain this apparent coincidence. 
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FIG. 3. The time dependence of the density p{t) for 
B-model with g = 3 and (from top) T = 2.9, 2.8, 2.7(critical), 
2.6, 2.5, 2.4, and 2.3 (L = 20000). Each line is an average of 
100 independent runs which starts from random initial con- 
figurations. Straight dotted lines have slopes corresponding 
to 5 = 0.286 and 0.5. 
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FIG. 5. The time dependence of the number of ac- 
tive sites N{t) for B-model with g = 3 and (from top) 
r = 2.9, 2.8, 2.7(critical), 2.6, 2.4, 2.2, and 2.0 (L = 50000). 
Each line is an average of about 10^ independent runs. The 
straight dotted line has a slope corresponding to 5 = 0.5. 
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FIG. 4. 
ability P{t) 
T = 2.9, 2.8, 
Each line is 
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The time dependence of the survival prob- 
for B-model with q = 3 and (from top) 
2.7(critical), 2.6, 2.4, 2.2, and 2.0 (L = 50000). 
an average of about 10^ independent runs. 
Straight dotted lines have slopes corresponding to 5 — 0.286 
and 0.5 respectively. 



The behaviour of p(t) (Fig. ||) is typical for PC univer- 
sality class. In the low temperature phase (T < 2.7) we 
observe the power law decay p{t) ~ t~^^'^ while at criti- 
cality (T = 2.7) p{t) - f-^ and we estimate S = 0.29(2). 
The dynamical Monte Carlo method confirms the PC 
criticality of this model |]l^. Indeed, P{t) scales with 
the exponent S' = 0.29(2) (Fig |) and at criticaHty N{t) 
remains virtually constant (Fig.^, which is in agreement 
with the PC value rj = 0.0. In our opinion, the above re- 
sults clearly indicate the PC criticality of B-model. This 
behavior can be qualitatively explained by considering 
the relation between our models and certain multispecies 



The relation with such models is based on the observa- 
tion that each domain wall can be at least approximately 
identified as a branching and annihilating random walker. 
Although precise mappings between absorbing-states and 
BARW models are usually complicated, one can argue 0| 
that most relevant processes are only long-lived ones 
which considerably simplifies the resulting BARW model. 
In particular Hinrichscn argued that dynamics of mod- 
els with two absorbing states should be approximately 
described by the following B ARW2 model : 

2X^0, X ^3X (2) 

This model is known to belong to the PC universality 
class Since there are several types of domain walls, 
mappings for models with more than two absorbing states 
require several types of random walkers. Hooyberghs et 
al. argued that typically in the corresponding g-BARW2 
models the following reactions should be included M : 



X 



Y + Z, 



Y + Z ^ X. 



(3) 
(4) 



The process (^) represents the formation of a domain of 
e.g., type C between domains of type A and B. The pro- 
cess (^) is its reverse and corresponds to the disappear- 
ance of the domain C. The g-species BARW2 model (p|- 
(^) differs from the one studied by Cardy and Tauber . 
However, the differences do not affect the critical be- 
haviour which is the same for both models. We expect 
that our A-model for q > 2 at the coarse-grained level is 
also described by q-BARW model with reactions (|^)-(Q). 

Let us now notice that in our B-model, the formation 
of the intruding domain C between domains A and B 
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is forbidden, hence processes of the type (g) are forbid- 
den. The g-species BARW2 model which corresponds 
to B- model is thus described by reactions and (||). 
The main point of our argument to explain the PC crit- 
icality of the B-model is the following: The suppress- 
ing of the processes (||) implies separation of time scales 
of parity conserving processes (|2|), which happen at the 
much shorter time scale than parity-nonconserving pro- 
cesses (^. As a result the g-BARW2 model spatially 
decomposes into single-species 'clouds'. The dynamics 
within each 'cloud' corresponds exactly to the dynamics 
of the q — 2 B-model (which is exactly the same as as 
for A- model), which at the coarse-grained level is given 
only by parity-conserving processes . Non-conserving 
processes (||) operate only when different 'clouds' col- 
lide [ po[ . Upon approaching the transition point do- 
mains coarsen and the distances between 'clouds' in- 
crease. As a result, non-conserving processes happen 
only on a very-long time scale. (Snapshots of Monte 
Carlo simulations qualitatively confirm formation of such 
single-species 'clouds'). On the other hand, the order 
parameter of the system, i.e., the number of particles 
is mainly determined by the parity-conserving processes 
(these processes determine the concentration of particles 
within clouds). As a result, the most relevant dynamics 
of the model is dominated by the parity-conserving pro- 
cesses which implies PC criticality. As we already men- 
tioned, these processes correspond to the coarse-grained 
dynamics of the q = 2 model which thus explains why 
transition temperatures for q = 2 and 3 are the same. It 
would be interesting to confirm these qualitative consid- 
erations with more sound theoretical arguments. 



III. CONCLUSIONS 



In conclusion, we have shown that the critical behavior 
of a model with q absorbing states is not only character- 
ized by the number of absorbing states but that details of 
the dynamics are also important. At the coarse-grained 
level our model is equivalent to a certain (7-BARW2 
model. It would be interesting to check whether sup- 
presion of processes (||) in q-BARW2 model leads to a 
similar change of the critical behaviour. 
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